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We investigate a class of microscopic systems in interaction with a macroscopic system 
in thermal equilibrium, following the construction of Dalibard, Dupont-Roc and Cohcn- 
Tannoudji (DDC). By considering self-adjoint operators as elements of Schwinger's Mea- 
surement Algebra (SMA), we construct statistical mean values of the relevant observables 
as matrix elements in a suitable operator basis, which correspond to the vacuum states of 
ThermoField Dynamics (TFD). 

1. Schwinger Operadors and TFD 

In SMA [1] an operator is denned as 

X = ^(a' | X | a")M(a',a"), (1) 

a', a" 

where (a' | X \ a") represents its matrix element in the a representation: 

(a' | X | a") = TrXM(a", a'). (2) 

The product of two general measurement symbols, referring to distinct sets 
of compatible observables A, B, C and D, satisfies the following composition 
rule: 

M{a', b')M{c', d') = (b'\c')M(c', d'), (3) 

where the number (b'\c') defines the statistical relation between the corre- 
sponding sets. 

The expectation value of a given property A in the b' basis is 

(A) v = a 'p( a 'i h ') = TrAM(b') = (b' | A | b') , (4) 

a' 

where we introduce the operator 

A = J2 a'M(a') = ^ a'5(a', a")M(a', a") , (5) 

a' a', a" 
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in the basis M(a') = M(a',a'). This implies that (a' | A | a") = a'5(a',a"). 
Now we define the statistical mean 



(A) = J2n(b')(A) b > = J>' | P A | a') = TripA), 

V a' 

with n(fc') > and n(fc') = 1, where 

p = W)M(b') = WW, b")M(b', b") . 

V b',b" 



(6) 



(7) 



Hence (V \ p \ b") = U(b')5(b', b") , with Trp = J2vW) = L In ther- 
mal equilibrium, by considering the basis M(n, m) associated to the number 
operator N in which its matrix elements are diagonal, 



we have 



where 



or 



N = ^2(n\ N \m)5 nm M(n,m) 

n,m 

= | N | n)M(n) , 

n 

p = J2n(E n )M(n), 

n 

U(E n ) = Z-V^" , 
p = Y,n(En)Sn m M{n,m). 



(8) 

(9) 

(10) 
(11) 



Introducing an auxiliary operator basis where a new fictitious operator N, 
corresponding to N, is diagonal, we can write 5 nm M(n, m) as 



so that 



(n | fh)M(n, m) = M(n, n)M(m, n) , 
V / n(K)n( J B m )M(n, n)M(m, n) 

,m 

£ V / n(K)M(n, n) £ \/n(£ m )M(m, m) 

n J L "i 

W))W)\ ■ 



(12) 



(13) 
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We may identify the measurement symbols inside the square brackets as 
composite states of the thermal vacuum of TFD [2]: 

| 0(/3)> ^ v 7 !^ \n,n), (14) 

n 

so that p acquires the status of a projection operator, as well as M(n,n) = 
\n)(n\ and M(rh,m) = \rh)(m\. Therefore, for a given observable F, 

Tr(pF) = {0(J3)\F\0(/3))=Tr(F\0(J3))W)\) = 

= Tr(^2y/U(E n )U(E m ) \m)(m\F\n)(r 

\ n,m 

= V^(E n )U(E m ) (n\m)(m\F\n) 

n,m 

= ^n(£„)HF|n>. (15) 



2. Radiation Considered as a Reservoir 

Let us consider the problem of a small system A interacting with a large 
reservoir 7Z. Following DDC construction [3], let 

H = H A + H R + V, (16) 

be the Hamiltonian of global system A + 7Z, where H A and H R are, respec- 
tively, the Hamiltonians of A and 7Z, whilst V stands for the interaction 
between them. In the interaction picture, the density operator of the global 
system obeys the evolution equation, 

from which we obtain 

ft+Ai 



p(t + At) = P(t) + - J [V(t'),p(t)]dt> + 

+ \Jh) I dt 'l ^(0.^(0^(0]]^- (is) 



By taking the trace with respect to TZ, we arrive at 

ft+At 



Aa(t) = a(t + At) - a(t) = — Tr R [V(t'), p(t)}dt' + 

inJt 

dt' jf Tr R [V(t'),[V(np(t")W, (19) 



ih 



I \ 2 r t+At 



where 



a(t)=Tr R p(t). 



(20) 



The interaction V between AeTZ will be taken as a product of an observable 
A of A and an observable R of 7Z: 

V = -AR. (21) 

Since the average value in a R of the coupling V(t) is zero, the leading con- 
tribution in f)19p stems from the two-time average 

g(t',t") = Tr R [a R R(t')R(t'% (22) 

If V is sufficiently small, and At sufficiently short compared with the evolu- 
tion time t R of cr, p(t) can be written in the form 

p(t)=Tr R p(t)®Tr A p(t), (23) 

where the contributions of the correlations between A and 1Z in the time t 
were neglected. The general idea is that the initial correlations between A 
and 1Z at time t disappear after a collision time r c <C t R . Thus, there exist 
two very different time scales, such that 

r c < At < t R , (24) 

bringing us to the 'coarse-grained' rate of variation for the system A. 

3. Master Equation for a Damped Harmonic Oscillator 

3.1. The Physical System 

We are interested in the case where the small system A is a one-dimensional 
harmonic oscillator of frequency ujq whose Hamiltonian is 

H A = fku (b j b+^), (25) 

where W and b are the rising and lowering operators of this oscillator. The 
reservoir 1Z consists of an infinite number of one-dimensional harmonic oscil- 
lators, of frequency uii, with ladder operators a\ and a^, so that the Hamil- 
tonian H R for 1Z is written as 

H R = Yl nw Ma + l)- (26) 
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We take a sesquilinear interaction between A and TZ of the form 

V = V^ = J2(v^a l + r ] *bal), (27) 

i 

where rji is the coupling constant between A and the i-th oscillator of TZ. 
3.2. The Master Equation 

We write the coarse-grained rate of variation, 

^ = -1± jf 00 d(t> - if') J*** dt'Tr R [V(t% [V(t'% 'A® ® a R ]], (28) 

with V given by (1271) . Changing to the Schrodinger representation, we obtain 
the following operator form for the master equation: 

%= ~ l[^b] + -r'[a^b} + -V'a 

- + /\)[tfb,a] + Tbatf + T'(tfab + bab ] ) . (29) 

In this equation [ ,] + represents the anticommutator and we have made the 
following definitions 

r = ~j[ ^2 I ^ I 2 - Ui) , r' = ^2 I ^ I 2 ( n ^(^o - Wt), (30) 

i i 

and 

/ia = v.p. (y ) , fua'=v.p.(y [ Vi |2( r } V (si) 

Here (rij) is the average number of excitation quanta of the oscillator i. If this 
number depends only on the energy of this oscillator, due to delta function 
in the second equation in (1501) . we have 

T = (n(ou ))T, (32) 

where (ti(ujq)) is the average number of quanta in the reservoir oscillators, 
having the same frequency ujq of oscillator A. If, moreover, TZ is in thermo- 
dynamic equilibrium, (n(ujo)) is equal to [exp fouio/hsT — l] -1 . 
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Alternatively, temperature effects can be incorporated from the very be- 
ginning, replacing the trace in (|28p by the thermal expectation value (|15|) . 
leading to 

^ = -~ 1°° d(f - 1") £ +At dt'wwit) ,[v(n <r A (t)]]W)), 

(33) 

and, consequently, to the master equation (1291) . In this case, the average 
number (rij) in (13"0"1) and (l3Tj) is given by a mean value in the TFD vacuum 
state, defined in (1T4")) . 

(n i ) = (0(p)\4aiW))> ( 34 ) 
after working out the thermal Green function 

g(r) = Tr[R(t)R(t-r)W))(0(f3)\) 

= (0((3)\R(t)R(t-r)\0((3)), (35) 

instead of the two-time average (f22|) . with r = t' — t". 

As a straightforward application of the formalism just presented, equation 
(f3"3"|) can be employed to derive the corresponding master equation (T2!l . in 
order to study a brownian particle interacting with a macroscopic system in 
the framework of equilibrium TFD [4] . 
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